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We have constructed the geometric phases emerging from the non trivial topology of a space- 
dependent magnetic field B(r), interacting with the spin magnetic moment of a neutral particle. 
Our basic tool, adapted from a previous work on Berry's phases, is the space-dependent unitary 
transformation U{r) which leads to the identity, W(r)^ S • B(r) W(r) = |B(r)| Sz, at each point r. In 
the "rotated" Hamiltonian "H, ^ is replaced by the non-Abelian covariant derivative — ^^(r) 
where A{r) = ihUK-§^U can be written as Ai (r)^^; -I- A2{r)Sy + A.-i{r)Sz- The Abelian differentials 
Afc(r) • dr are given in terms of the Euler angles defining the orientation of B(r). The non-Abelian 
field .4(r) transforms as a Yang-Mills field, however its vanishing "curvature" reveals its purely 
geometric character. We have defined a perturbation scheme based upon the assumption that in 
H the longitudinal field A3(r) dominates the transverse field Ai.2(r) contributions, evaluated to 
second-order. The geometry embedded in both the vector field A3(r) and the geometric magnetic 
field B3(r) = ^ A A3(r) is described by their associated Aharonov-Bohm phase. As an illustration 
we study the physics of cold ^^^Yb atoms dressed by overlaying two circularly-polarized stationary 
waves with orthogonal directions, which form a 2D square optical lattice. The frequency is tuned 
midway between the two hyperfine levels of the (6s6p)^Pi states to protect the optical B(r) field 
generated by the lattice from the dressed atom instability. The geometric field B3(r) is computed 
analytically in terms of the Euler angles. The magnitude of the second-order corrections due to the 
transverse fields can be reduced to the percent level by a choice of light intensity which keeps the 
dressed atom loss rate < 5 s~^. A second optical lattice can be designed to confine the atoms inside 
2D domains where B3(r) ■ i > 0. We extend our analysis to the case of a triangular lattice. 



PACS numbers: 03.65.Vf,37.10.Vz 

I. INTRODUCTION 

In this paper we shall be interested in the spatial geom- 
etry associated with the time- independent Hamiltonian 

describing the quantum evolution of a non-relativistic 
neutral particle of spin S and magnetic moment 75 S, 
interacting with a scalar potential V{v) and a static non- 
uniform magnetic field B(r). 

In a recent work concerning Berry's phases Q gen- 
erated by arbitrary spins non-linearly coupled to time- 
dependent external fields 0], we have found that the 
geometry of the Hamiltonian parameter space is more 
clearly exhibited if one uses a rotating frame method 
instead of the standard approach involving space-time 
wave functions within the adiabatic approximation. The 
Coriolis effect generates in the rotating frame a time- 
dependent linear spin coupling AH{t) — — 75AB(t).S. 
The Berry phase is generated by the longitudinal effective 
field ABz{t) and the non-adiabatic corrections, governed 
by the transverse fields ABx^y{t), are readily obtained by 
a second-order perturbation calculation. 



The purpose of this paper is to use a similar approach 
to construct the geometrical phase associated with the 
non-trivial spatial topology of the magnetic field B(r). 
We are going to apply to "H a local unitary transfor- 
mation, which makes diagonal the spin coupling at each 
point r. The transformed kinetic term is obtained by 
replacing the gradient ^ by the non-Abelian covariant 
derivative, D = |: - j-Ak{r)S^ (with = S*^, = 
Sy, — Sz), in which is encapsulated all the geome- 
try of B(r). By analogy with Berry's phase it is possible 
to assume, within well-defined conditions, that the longi- 
tudinal field A3(r) gives the dominant contribution, the 
transverse ones, Ai_2(r), leading to second-order correc- 
tions. This is a kind of Born-Oppenheimer approxima- 
tion, playing, here, the role of the adiabatic approxima- 
tion. The role of the Berry phase is, then, taken up by 
the Aharonov-Bohm phase associated with A3(r). 

A simple example of this problem, where the spatial de- 
pendence of B(r) results from a constant field gradient, 
was studied experimentally by W. Phillips and coworkers 
In the case where B(r) is a periodic field, Ti. could 
be used, for instance as a model to describe the scatter- 
ing of thermal neutrons by magnetic materials [3| . There 
has been recently a regain of interest in this model for 
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describingthe evolution of cold atoms within an optical 
lattice UHl- In this case, B(r) stands for the "effective" 
magnetic field generated by an ac-Stark effect involving 
the electric field E(r) associated with the coherent laser 
field used to construct the optical lattice. Two realis- 
tic and simple examples involving S = 1/2 cold atoms 
trapped in a two-dimensional (2D) optical lattice will be 
discussed in the last section of this paper. 



II. NON-ABELIAN GAUGE FIELDS 
GENERATED BY THE B(r) GEOMETRY 

A. The " rotated frame" approach 

We are going to extend the method we have used in 
the case of time-dependent fields to the present context, 
by introducing at each space point r a "rotated frame" . 
The change of frame is defined by writing the vector state 
as 4'(r) = U{v)'^[r). The space-dependent unitary trans- 
formation Z//(r) is designed to line up, at each space point, 
the magnetic field B(r) along the spin quantization axis: 

U\v)S-B{v)U{v) = \B{v)\S,. (2) 

Let us ignore, for a moment, that hl{v) does not com- 
mute with the kinetic energy term. For a periodic mag- 
netic field, one would then have to solve, for each spin 
component m, a wave equation involving a scalar poten- 
tial, like in the historical experiment of Stern and Ger- 
lach. One would then come down to solving standard 
problems of solid state theory. 

To proceed, let us introduce the rotation, 7?.(0, 8) ~ 
R{z,(p{r)) ■ R{y,9{r)), where R{n,x) is the rotation of 
angle x about the unit vector h and <j){r), 9{r) the two 
Euler angles which define the direction of the B(r) field. 
More precisely, TZ{(j){r), d{r)) is chosen in order to satisfy 
the following equation, 

B(r) = |B(r)|7e(0(r),0(r)).z. (3) 

We should keep in mind that there is not a unique way 
to bring z along the direction of B(r). This property 
can be associated with the gauge invariance of the non- 
Abelian gauge fields that we are going to introduce. It 
is easily verified that U{r) given below coincides with 
the unitary transformation associated with the rotation 
7^(0(r),0(r)): 

W(r) = cxp(^ ^(r)). exp(^ Sy 0(r)). (4) 

Using group theory arguments, one can derive the basic 
relation (see ^]), 

Wt(r)SW(r)=7^(0(r),0(r))•S. (5) 

This leads to the set of identities: U''{r)S • BZ^(r) = 
(7^(0, 0)S)-B = S- (7^-l(0, 9) B) = |B(r)| S,. 



B. The geometric non-Abelian gauge fields 

Our next step is to write within the rotated frame the 
eigenvalue equation, "H \E'„ = now taking into ac- 

count that U(r) does not commute with the kinetic en- 
ergy operator ym' making simple manipulations, we 
arrive at the eigenvalue equation written in the rotated 
frame: ■H^„(r) = i?„^'„(r), where the rotated Hamilto- 
nian H = U\r) % U{v) is given by: 

n = ^p2 + nr)-75|B(r)|5„ (6) 
P - -in^-inu\v).{^U{v)). (7) 

In reference (Q sec. 4.2), we have made a similar calcula- 
tion in the case where the Euler angles are time- instead 
of space-dependent. The needed result is obtained by do- 
ing the replacement: >0. This 
leads to the following expression for p: 

d 

Ai(r) = -sine'(r) ^0(r), 

A2(r) = — 0(r), A3(r) =cos%)— 0(r). (8) 
or or 

Note that A3(r) 5*2, etc. are tensor products of a vector 
operator in the geometric space by a vector operator in 
the spin space. A remarkable feature of the above results 
is also the fact that the three fields Ai(r) do not depend 
on the value of = S{S + 1)%^. This follows from 
the fact that the evaluation of p relies only upon the 
SU{2) Lie algebra: [S\ S^] = ihtijk . A simple direct 
evaluation of the above formulas has been performed for 
S" = 1/2 using the fact that U{r) is given by a 2 x 2 
matrix, linear with respect to S. From now on, we shall 
often use the notation = S^, = Sy, = Sz 
in order to stress that the S'' will be treated as scalar 
objects under ordinary space rotation. It is of interest to 
incorporate the three fields Aj(r) into a "non-Abelian" 
classical field A{r), a well-known concept used in other 
fields of physics: 

^(r) = Afe(r)5'= = ihU\r).i^U{r)). (9) 

The Hamiltonians considered in this work have not any 
simple invariance properties under rotation involving the 
total angular momentum of our neutral particle J = 
L -t- S. Using a fashionable expression, our spin will be 
treated as a "colored" spin. 
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C. Non-Abelian covariant derivatives and gauge 
invariance 

In analogy with what is done for Abehan fields, we de- 
fine the basic concept of covariant "derivative" as follows: 

We are now going to require that D.\I'(r) transforms 
under any space-dependent SU (2) unitary transforma- 
tion exp(^Afe(r) S^) in the same way as the wave func- 
tion $(r), i.e. $(r) ^ $A(r) = exp(iAfc(r)5''=) $(r)). 
The above requirement will be satisfied provided that 
the non-Abelian field ^(r) is subjected to a well-defined 
non-Abelian gauge transformation: ^(r) ^A(r). The 
covariance condition means that the operator D should 
transform locally as a vector operator under any space- 
dependent SU (2) unitary transformation: 

|: - J^A = exp( jAfc5'=) D exp(- jAfc^'^). (11) 

The non-Abelian gauge transformation A{y) — > .4a (r) 
is then easily derived by introducing in the above formula 
the expression of D given in Eq. (flOt , 

^A(r) = exp(jAfc5'=)(^exp(-jAfe5'=)) 

+ exp(jAfe^'=)^(r) exp(-^Afc5'=). (12) 

It is found to agree with the formulas given by S. Wein- 
berg for general non-Abelian fields. If the space- 
dependent matrix Ak{r)S'' is replaced by an ordinary 
scalar function of r, one recovers the usual Abelian gauge 
transformation. 

An important physical concept for non-Abelian gauge 
fields is the "curvature" tensor, involving covariant 
derivative commutators, Tij{r) = i[Di,Dj], or 

d d 

A crucial point, here, is that the non-Abelian field A{r) 
given by equation ([T^ can be expressed under two 
forms, which are equivalent as a result of Li(v) unitar- 
ity: ^(r) = ^hUHr) . {-§-U (r)) = -th{^W {v)) .U{v)). 
Using this identity to rewrite -Fijir) one arrives, after 
some algebraic manipulations, at the remarkable result: 
Tij{r) — 0. This shows, clearly, that .4(r) is not a self 
interacting Yang-Mills field like the fields used in Particle 
Physics, but a purely geometrical field associated with a 
local change of space coordinates. There is a similarity 
with the situation encountered in the tests of the "equiv- 
alence principle" which involve "flat" space-time metrics. 

A possible application of the non-Abelian gauge trans- 
formation is the choice of the Landau gauge for AA3(r) 
by writing A/^s^ = x ■ AA3(r) = without affecting the 
magnetic coupling. It is easily seen from Eq. (11), that 



the wanted gauge transformation has to satisfy the fol- 
lowing constraints: Ai = A2 = , ^ A3(r) = —Asxi'r)- 
An explicit expression A3 (r) is obtained by a quadrature 
over y : A3(r) = — A^^ivj y, z) dv. With this choice of 
A the Hamiltonian H is obtained from H by performing 
the replacement Ai(r) — > Ai(r), 

Ai(r) = cos(A;3)Ai(r)-f sin(A3)A2(r), 
A2(r) = -sin(A3)Ai(r) +cos(A3)A2(r), 

A3(r) = {Asy{r) + -^As{r))y. (13) 



III. THE AHARONOV-BOHM PHASE FROM 
THE "LONGITUDINAL" GAUGE FIELD 

Although the above developments are valid for any spin 
values, from now on, for sake of simplicity, we shall limit 
ourselves to the case of a S — ^ particle in a 2D space. 
So far we have made no approximation. In analogy with 
the method we have used for Berry's phases generated by 
non-linear spin Hamiltonians [9], we are going in sec. IV 
to construct a perturbation expansion, using as starting 
point the diagonal part Ho of the Hamiltonian T-L, 

{...} denoting the anticommutator between operators. 
The above approximation is often referred to as a Born- 
Oppenheimer (BO) approximation. It supposes that the 
evolution of the spin wave function is slow as compared to 
that of the spatial wave function describing the particle 
internal state (just like in molecular physics, the motion 
of the nuclei of a molecule is generally slow in comparison 
with that of the atomic electrons). The validity of this 
approximation in the present context is dicussed later 
on (see sec. IV). This is one advantage of the present 
"rotated" frame approach, similar to our treatment of 
the adiabatic approximation j9i], to provide a systematic 
perturbation method to calculate the corrections. The 
left over perturbation Hamiltonian is given by the non- 
diagonal contribution of 

^1 = • (Ai(r)5, + A2{v)Sy) + h.c. (15) 

Let us discuss the physics content of Eq. [14] for a given 
eigenvalue of Sz/h, to = ±i. The three first terms de- 
scribe a spinless particle moving in a scalar potential. 
The terms of the second line represent the coupling of this 
particle of effective electric charge to, interacting with the 
Coulomb-like potential — 7s|B(r)| and a magnetic vector 
potential given by A3(r). The diamagnetic term A3(r)^ 
is incorporated into the third term of the first line. 
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FIG. 1: The geometric magnetic field Bs^ and the confinement potential inside one elementary cell of the optical lattice 
represented in the "rotated frame" . The Left graphic gives a perspective 3D plot of Be™, (r) , the electromagnetic field equivalent 
of the geometric field B3(r) = ^ A A3(r), as a function of fca; and ky ik being the wave number of the optical "dressing" field 
generating the non-Abelian gauge fields). The geometric magnetic field -B32 exhibits two positive maxima within the plaquette 
defined by 7r/2 < kx < Stt/2, —n/2 < ky < tt/2. A symmetric pair of negative minima Bsz appears within the plaquette 
obtained by performing the translation a: ^ a; + tt, as predicted by the formula (|27|) . In the Right graphic we have plotted Bsz 
together with the total potential Vc{r) — V(r) + VV(r) — 7/ h/2 |B(r)| acting upon a state with ~ fi/2, drawn in dark red 
using as unit 7/ 7i Bo- ( For the sake of visibility Vc(r) has been shifted up by one unit.) As expected, the added potential VV(r) 
does prevent the cold atoms from exploring the regions where the sign of Bsz is negative. Indeed, the pair of negative minima 
Bsz of the Left are clearly not accessible for atoms below a well defined temperature. It is of interest to note that the two 
graphics are invariant under a rotation of vr around a vertical axis centered at the point kxo = Sn/2, kyo = vr/2, leading to the 
coordinate transformation, x^x' = 3n — x,y^y' — n — y which leaves invariant both -B3z(r) and VV(r). The values of these 
quantities for the full 2D space are obtained by performing the set of translations: x—^x' = 2ni-K + x,y-^y' — 2n2Tr + y 
where ni and n2 are arbitrary integers. 



Ignoring for a moment Hi, we would like to con- 
sider a two-slit Aharonov-Bohm interferometer experi- 
ment involving two beams described by wave packets, 
solutions of the Schrodinger equation associated with "Ho- 
lt is convenient to describe the two interfering beams by 
two Feynman path integrals associated with the classical 
paths Ci and C2 enclosing a 2D finite region TZb where 
the "geometric magnetic" field is confined: 

B3(r) = |-AA3(r)=sin(^?(r))(^0(r)^0(r)- 

^^(r)|:.(r)).. (16) 

Following Feynman [ll|, the phase difference 

^ABim) = -m(^<j> A3(ri)dri-^ A3(r2)dr2), 

(17) 



is factored out of the functional integration and is given 
explicitly by the 2D integral: 

$as("i)=to/ dxdy B3^{x,y), (18) 

where mfi plays the role of the electric charge. The 
Aharonov-Bohm phase [13] $ab(w), can be viewed as 
a spatial extension of Berry's phase associated with the 
non-trivial geometry of the physical B(r) field, which is 
described by the geometric magnetic field, B3(r). 

The above calculation has been performed in the ro- 
tated frame. To return to the laboratory frame, we shall 
use the fact that, in a typical Aharonov-Bohm experi- 
ment, the paths Ci and C2 are drawn in regions where 
B3(r) = 0. Assuming that the Feynman phase factors 
A3{ri) dvi have been factored out, the wave func- 
tions relative to the paths Ci can be written in the ro- 
tated frame as = '4'i,m{i^, t) <S) Xm- The spatial wave 
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functions '(pi^mi^^t) obey the Shrodinger equation associ- 
ated with the Hamiltonian ■Ho(A3 =Q,Sz = mTi) given 
by Eq. (jl4p : Xm is the S — 1/2 spinor wave function, 
('5m.i/2i '^m,-i/2)- When going back to the laboratory 
frame the spatial wave function ipi^rn{i^Tt) is unchanged, 
whatever r, t, since U (r, t) acts only on the spin, while the 
spinor, Xi,m(r) = U{r,t). Xi,m has been rotated. How- 
ever, at the interference point the spinor wave function 
interference in the laboratory and in the rotated frame 
is identical, (x2,m|xi,m) = {X2,ni\l^H^)-^i^)\xi,m} , as a 
consequence of the unitarity relation U^{r)M{r) = 1. In 
more physical terms the spins rotate differently along the 
two paths Ci and C2 , but their scalar product is preserved 
at the interference point. In conclusion: $^^(771) = 
$As(m). 

In order to clarify the relation of ^AB{m) with the 
standard Berry's phase, it is instructive to discuss briefly 
the case of a space and time-dependent B(r, t) field. 
As before, let us construct the wave function in the 
frame rotating in space-time: ^(i) — {r,t)'i>{t), 
where ^'(t) is a solution of the Schrodinger equation 
associated with the time-dependent Hamiltonian Hit) 
obtained by giving a space- and time-dependence to 
the magnetic field B. The unitary transformation 
Zi(r,i) = exp(^5',0(r,t)).exp(^S'j,6'(r,t)) brings the 
B(r,t) field along the z axis. The time evolution of the 
"rotated" wave function is then governed by the Hamil- 
tonian n{t) = ^{-ifif^ - Au{v,t)S^f - Vu{v,t)S'' + 
Vir) - 7s|B(r,t)|S'^. The scalar fields Vk{v,t) are ob- 
tained from the vector fields (r, t) by replacing in 
Eq. © ^ by For the moment, let us remove the r 
dependence and ignore the transverse potentials Vi,2(i), 
governing the non-adiabatic corrections 9]. Berry's 
phase is closely related to the dynamical phase contri- 
bution generated by —Vz{t)Sz in the rotating frame @: 
<i>_D3 = /g dt m cos 9 {t) 4>{t). By returning to the lab- 
oratory frame, ^(i) ^'(i), one recovers the standard 
Berry's phase formula: fi{m) = dt m cos 6 {t) — l)<j){t). 

The creation of effective magnetic fields in 2D optical 
lattices was suggested by Jaksch and Zoller [l3|. They 
deal with a pseudo spin i associated with two internal 
atomic states. The effective fields are resulting from the 
spatial dependence of the Rabi frequencies induced by 
Raman lasers designed to transfer an atom in a given 
internal state to a nearest-neighbor lattice site associ- 
ated with a flipped internal state. In addition, there is in 
the Hamiltonian a dynamical contribution which specifies 
the direction of the atomic jump. This requires, in prac- 
tice, additional space-dependent potentials [14]. Here, 
we have adopted a totally different point of view: we 
have focused upon continuous space-dependent "dressed" 
Hamiltonians, having a linear spin coupling. We have de- 
veloped a precise formalism showing how the non-trivial 
spatial topology of our spin Hamiltonian implies the ex- 
istence of an Aharonov-Bohm phase, generated by a ge- 
ometric magnetic field B3(r). The role of the electric 
charge is taken up by the eigenvalue of Sz/^, m = ±i. 



IV. THE TRANSVERSE GAUGE FIELDS 
CORRECTIONS 

We are going to present a brief analysis of the effects 
or the transverse gauge fields Ai^2(r) upon the physics 
associated with the approximation 'H « 'Ho, taking as an 
example the calculation of the Aharonov-Bohm phase. 
We shall use a perturbation approach using as starting 
point the eigenstates of Hq associated with a given eigen- 
value of to: Ho'^nmir) = Enm'^nmir). The eigenenergy 
corrections appear to be of second-order with respect 

to Hi: A2En„^ = l(*n,m|^l|* 

,m 

En^^m)- Assuming that the magnetic coupling energy, 
~7S'S'z|B(r)|, is dominant in "Ho allows us to write: 
En,m — En,^ra as thc quantum average of —2m^s |B(r)|. 
The energy shift /S.2Enm is then approximated by the 
rather simple expression: 



2TOn7s(*„„|(|B(r)|)|^', 



(19) 



Let us now turn to the "Hi corrections to the geomet- 
ric phase $^b(to). In a typical atomic interferometry 
experiment, the convergence of the two paths Ci is im- 
plemented with a localized "mirror" device. If the gauge 
vector field line integral is factored out, the wave pack- 
ets $i(r, associated with d satisfy, outside the "re- 
flexion" regions, the Schrodinger equation governed by 
"Ho (A3 = 0), a real operator having real eigenfunctions 
'^n,m with eigenvalues £n,m- Let us concentrate upon 
wave packets relative to the converging sections of the 
paths, $,= (r, t) with < t < T, where the times -T, 0, T 
are respectively the emission, reflexion and interference 
times. At these three times occur transformations de- 
scribed within the sudden approximation. Using stan- 
dard time-dependent perturbation formalism, the first- 
order T-Li corrections to the converging wave packets 
reads: 



m{t) 



2 z ^($fi,_„,|Hi|$„,„0 ($„,™|$,"(0)) 

n.n 

Sin((£:„,„ - £.n,-ni)t/h) 



Since only the wave packets with identical m can in- 
terfere, the lowest-order interference correction reads: 
( S^liT) I 5n{T) ) = ( 5^1{T) I 5n{T) ) cM^^^ab)^ 
where A^\b phase difference accumulated along 

the converging paths. A$^^ — rn(^f^^ A3(ri)(iri — 

(Ci — > C2)). (Note the change of sign with respect 
to ^AB-) To proceed we make two approximations. 
As before in the evaluation of A2i?nm we shall write: 

(Sn^rn - £n,-m) = 2 TO 75 ( ™ | ( |B (r) | ) | ^) . Then, 

we shall drop the oscillating sin products and replace the 
sin^ by their average 1/2. This leads to the following 
expressions for 2 = (^*2(^)l^^i(^)) ^^d 77„,,„, 

'^1^2 « ^?7n,m ($2(0)l<I'„.™>('fn.™l*^(0)), 
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(7s?i(*nm|(|B(r)|)|*„„0)2^ 



(20) 



which are closely related in magnitude to the energy shift 
^2Enm, calculated previously (Eq. [T9| . 

It is of interest to give an explicit expression of 
($„,m|H?|$„,™), by writing Hi = ELi ^i,fe ^''^ where 
ili^k = {p , Afc(r)}/(2M). Using the S'' algebra for 
S = 1/2, one gets the more compact expression. 

Pulling out the Afc(r) field from the matrix element of 
"Hi.fc one gets a formula useful for an explicit evaluation: 



E 



Z'^''"^*"'" ^k,^Ak^j^ (21) 



Using the above results which make no assumption 
about the origin of B(r), it is possible to give an order 
of magnitude of rjn^m in the particular case of a periodic 
field B(r) associated with a single wave number k. Look- 
ing at Eq. ([5]) one sees clearly that the non-abelian fields 
Ai(r) scale as k and their divergence, ^ • Ai(r), as . 
In an explicit computation to be presented later on, it 
appears that the terms involving the field divergence is 
the dominant one. Assuming that the magnetic coupling 
is the dominant term in 'Ho, one arrives at the estimate. 



< 



1 / (nk)V(2A/) \2 ^^^^ 
wth ■nmax = -{- — r-T^-n — -] , 22 



which can also be viewed as a scaling law. In the case 
where the periodic field B(r) contains an effective opti- 
cal magnetic field Bopt(r) generated inside a 2D optical 
lattice, (?ik)^/ (2 AI) is the recoil energy of the cold atom 
following the absorption of one photon of momentum k. 



APPLICATION TO "^Yb COLD ATOMS IN A 
2D OPTICAL LATTICE 



We now want to illustrate the above theoretical analy- 
sis with a simple realistic example. We choose the case of 
cold neutral atoms with F = I = ^ angular momentum, 
"dressed" by a second-order Stark effect [l^ involving 
the wave electric field E(r) of an optical lattice. The light 
wave directions and polarizations are selected to generate 
an effective magnetic field Bopt(r) which has the lattice 
periodicity. 

Much more ambitious experimental programs involv- 
ing 2D optical lattices are aiming at a simulation of actual 
condensed matter problems [14lll6l - [l8| . like the electronic 
properties of graphene. 



E 



A. The dressed atom Hamiltonian T-L 

Let us assume that the "dressing" optical electric field 
E(r) is obtained by the coherent superposition of two sta- 
tionary waves, directed along x and y. Each one results 
from the interference of a circularly polarized laser beam 
(wave vector k, helicity ^) and of the beam refiected at 
normal incidence, all beams carrying along their propa- 
gation axis the same angular momentum per photon, 

{r)—£^—j={{l + z^z)(exp ikx — exp —ik x) 

=(i -f i^i)(expiA:y — exp— ifcj/)|. (23) 

This configuration has been selected in order to make 
the space dependence of iE* AE appearing in B(r) both 
simple and symmetric under the x ^ y exchange (see 
Eq. (|25l) below). In particular the relative phase of the 
two stationary waves plays an essential role. To avoid 
the singularities of the non-Abelian gauge fields directly 
associated with the geometry of the optical "effective" 
magnetic field Bopt, a uniform dc magnetic field Bq has 
to be applied, with the same x y symmetry as Bopt, 
Bq = (sinu {x + y)/y/2 + cosu z)Bo. 

We consider non-interacting ^^^Yb atoms, with nu- 
clear spin i, trapped in this 2D optical lattice and ig- 
nore many-body effects. The laser frequency is detuned 
from the transition at 555.8 nm between the ground state 
(4/",6s2)iS'o and the first (4/", 6s6p)3Pi excited state. 
In the case of a single transition between two states of 
angular momentum i , the optical field generates a space- 
dependent effective dc magnetic field Bopt which follows 
from the identity, valid for S — 1/2, 

(S • E*)(S • E) = — E* • E + i-E* AE- S. (24) 

The magnitude of Bopt scales as the atom-field coupling 
ri^/(5, where Tifl — 2£d is the Rabi frequency of the 
transition, d = {^Po, ni = 0\dz\^So) its electric dipole 
amplitude and 6 the laser frequency detuning. Actually, 
the ^^^Yb isotope has two hyperfine (hf) components, 
F = i and = | with relative strength i and | and 
hyperfine splitting Whf/n = 5.939 GHz It is impor- 
tant to adjust the laser frequency midway between the 
two hf lines, {i.e. opposite detunings, hS = ±Whf/2, 
for the two transitions). This particular choice has the 
great advantage of canceling out the imaginary part of 
the effective magnetic field associated with spontaneous 
emission (see Appendix). Only the number of dressed 
atoms, irrespective of their nuclear spin orientation, is 
subjected to a decrease. With this choice the contribu- 
tions to Bopt coming from each transition happen to be 
equal. 

It is convenient to write Bopt in a way which makes 
its coupling to the -"^^^Yb ground state exactly similar to 
the real magnetic field coupling. Since it originates from 
an ac Stark shift, it has, however, no reason to have a 
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FIG. 2: Illustration of the geometry embedded in the effective 
field B(r) by a perspective view representing the 2D surface 
S, immersed in a 3D Euclidian space. 5 is defined by the 
parametrization: xi = Bx{x,y)/Bo ,yi — By{x,y)/Bo , z\ — 
Bz{x,y)/ Bo- Performing the simple change of coordinates: 
X2 — xi — 1/2 — sm^{kx);y2 = yi — 1/2 — a'm^ (ky); Z2 — 
z\ — — — sin(fc a;) sin(fc y), one finds that S is the in- 

tersection of a cone of equation: 22 — X2 1/2 with the cubic 
volume defined by < 3:2, 3/2 < 1, together with the condi- 
tion — 1 < 22 < 1. The "rotated" frame method is associated 
with the mapping of S onto a line segment along the zi axis. 
The cone tip can be associated with the saddle point appear- 
ing in Fig. 1 between two bumps of the geometric field B3. 



vanishing divergence. From the calculation presented in 
the Appendix, Eqs. ([551 we obtain the expression, 



B 



8d2 



opt 



iE*(r) AE(r) 



37/ hWhf 

= Bq p sin^(fc x)x + a; — > y — sin(fc x) sin(fc y)z^ , 

_ 32 £^ _ 8 nn^ 

^~ iWhfhjiBo ^ SWhf-fiBo' ^ ' 

Then, the dressed atom Hamiltonian is expressed as: 

^ = ^-7/(Bo + B„pt(r)).I + V(r), 

V(r) = ~n-fiBo p{sm'^{kx)+sm^{ky))/'i, (26) 

where the scalar potential V(r) is associated with the 
scalar terms cx E* • E (see Eq. ([35)) V 



B. Explicit evaluation of the phases associated 
with the magnetic field geometry 

In the following, the AB phase will be considered as a 
diagnosis. It serves to test for the presence of a vector 



potential A3(r) able to generate the non-local quantum 
physical effects associated with the magnetic field B3(r). 

We have all we need to construct the non-Abelian 
gauge fields given by Eqs. [H] The angles ^(r) and 0(r) 
appearing in the relevant unitary transformation: U (r) — 
exp(— |tT2 (^(r)). exp(— I (Ty 6'(r)) are given by the orien- 
tation of the total magnetic field, B(r) = Bo(r)+Bopt(r), 



tan(/)(r) — 



By{r) 



cos0(r) = 



S.(r)' \B{r)l 

The graphs of Fig. [T] and [2] have been obtained for the 
typical case u = 7r/4, ^ = p = 1. The geometric mag- 
netic field B3(r) can, then, be obtained directly from the 
compact geometric formula (1161) . The algebraic structure 
of 0(r) and 9{r) suggests that B3(r) can be written under 
the symmetric form, 

B3(r) = — z fc^ cos(fc a;) cos(fc y) J^( sin(fca;), sin(fcy)), 

T{a,b) ^ {a^ + 2V2ab + b^)G{a,b)-^, 



G{a,b) 



1 



V2ab + b^ + b'^ + a^b'^. (27) 



By looking at Fig.l, it is clearly of interest to compute 
the phase $geom(l/2) as the flux of mB3(r) through the 
plaquette: 7r/2 < kx < 3tt/2, —tt/2 < ky < tt/2. Using 
the above explicit expression of B3(r), and performing a 
precise numerical quadrature, one obtains <I>geom(l/2) = 
0.499284. This pure number is a characteristic of the 
geometry associated with B(r) for the values of the geo- 
metric parameters u = 7r/4, ^ ^ p = 1. 

It is of interest to construct, from the geometric field 
B3(r), an electromagnetic field which would induce the 
same physical effect {e.g. the same Aharonov-Bohm 
phase for a particle with elementary electric charge e). 
This can be achieved by comparing the covariant deriva- 
tive of Ea. ([TU|) . in the hmit Ai(r) = A2(r) = (namely 
^ — imA3(r), mh being the eigenvalue ±i of Iz), 
with the corresponding expression in the electromag- 
netic case: ^ — i-| Ae,„(r). A simple inspection leads 
to: Aem(r) = TO|-A3(r). This implies Be„(r) = 
±i(3^)2|(B3(r)fc-2)^ where A(Yb) = 0.5558x10-^771 

is the wavelength of the ^5*0 — >■ '^Pi "'^''^Yb transition. Us- 
ing the values of e and h, one gets the electromagnetic 
equivalent in MKSA units of the geometric field B3(r), 
Be™(r) = ±0.0420287 x 10"* Tx(B3(r) fc-^). 

It is important to note that in the expression of B^z, 
Eq. ([?f)) . J^(sin(fca;) sin(/c7/)) is a positive definite func- 
tion of x,y. This suggests the possibility to confine the 
Yb atoms in 2D domains where B3(r) • z > 0, by adding 
to the H Hamiltonian an extra scalar potential. 



>V(r) — h-fj Bq p cos(fc x) cos{ky), 



(28) 



where p is a dimensionless positive parameter to be ad- 
justed. Let us write, in the "rotated" frame, the confining 
part of the Jio Hamiltonian associated with the dressed 
(4/" 652)15-0, 7=1 ground state of ^^^Yb, 



He 



P 

2M 



(V(r)-f W(r)~7, |B(r)|/,^ 



(29) 
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FIG. 3: The scalar potential Vc(r) (dark red) and the geomet- 
ric field (light blue), within one quarter of the elementary 
cell relative to the potential W(r). The geometric phase rel- 
ative to this plaquette cancels, but when W(r) is increased 
adiabatically, the trajectory of the wave packet is clearly mod- 
ified: the cold atom is attracted towards the B^z > regions 
and gets submitted to a Laplace-type force. The quantum 
circuit no longer satisfies the conditions required for the ob- 
servation the AB effect. Therefore the cancellation of the 
purely geometric phase no longer implies the vanishing of the 
magnetic flux through a physical circuit, hence cannot be con- 
sidered as a negative test of the Ba^ fleld relevance. 

We have omitted the gauge vector fields (r) contribu- 
tions, which play a limited role in the confinement. In 
the rhs of Fig. [1] we have plotted in red, using ^jBTi 
units, the total potential acting on the hyperfine sub- 
state Iz = ?i/2, 

K(r) - V(r) + W(r) - 7/ n/2 |B(r)|. (30) 

We have found that the choice /i = 2 leads to the wanted 
effect. On the same graph, we have displayed in blue 
twice the dimensionless quantity B3(r)fc^^. It is clear 
that, for an appropriate temperature, cold ^^^Yb atoms 
can be confined in connected domains of the 2D-space 
where the effective magnetic field satisfies B3(r) • z > 0. 
For this purpose, the experiment should make use of two 
optical lattices having the same periodicity, but each one 
fulfilling a different role. The green lattice (555.8 nm) 
provides the topology of the effective field which gener- 
ates the geometric non-Abelian field and the associated 
AB phase. The confinement role is attributed to the 
second lattice, spin-decoupled, hence not affected by the 
"rotated" frame transformation. This provides the pos- 
sibility to adjust the atom position according to their 
temperature. 

One way to generate the scalar potential yV(r) is to 
rely upon the scalar ac Stark effect with a large detuning. 
This may involve two stationary waves linearly polarized 



along z, having wave vectors ki = k{x + y)/2 and k2 = 
k{x - y)/2 with |ki| = jka] = fc/v^, i.e. X = 786.0 nm. 
We make the further assumption that there is no phase 
coherence between the two stationary waves, in other 
words, they do not interfere. With an appropriate phase 
and intensity matching of each stationary wave, the light 
induced potential, up to a constant factor, is given by: 
cos( k {x + y)/2) -|-cos( k {x — y)/2) = 2 cos{kx) cos(fcy). 
It has clearly the right space dependence for the wanted 
potential and the wanted sign is provided by the neg- 
ative red detuning. The light intensity can be used to 
adjust the potential depth without affecting the stability 
of the dressed atoms, thanks to the large magnitude of 
the detuning ~ lO^Fp. 

An important property of the fiux of B3 is its cancella- 
tion when it is taken through a full elementary cell. This 
is an obvious consequence of the periodicity of the gauge 
vector field A3 (r) generated by the periodic optical field 
Bopt(r). From a look at Fig.l, it is clear that the symme- 
try properties of B3 imply that its flux ^geom through 
the plaquette now associated with the half- elementary 
cell, defined by 7r/2 < kx < 57r/2 , -7r/2 < ky < 7r/2, 
vanishes. The vanishing of $geom - a purely geometrical 
result - is not affected by the scalar potential W(r) which 
has no effect upon the A3(r) gauge field. In the limit 
where the kinetic term dominates Vc(r) in the Hamilto- 
nian Jic, it is possible to build wave packets which propa- 
gate along the plaquette perimeter. The AB phase, ^ab, 
around the half-elementary cell will vanish since it coin- 
cides with ^geom- However, if we assume that >V(r) is 
increased adiabatically, the wave packets are pushed to- 
wards the regions where B^z > 0. Depending on their ki- 
netic energy, they may even be confined there. In these 
conditions, the phase accumulated along the deformed 
circuit will differ from $geom for two reasons: first, the 
B^z fiux is no longer guaranteed to vanish; second, the 
wave packets propagate in regions where {B^z \ > which 
makes the situation incompatible with physical require- 
ments for measuring the AB effect, see Fig. [3l There- 
fore the vanishing of the B^z flux through this particular 
purely geometric loop cannot be considered as an evi- 
dence for the irrelevance of the B3 field in presence of a 
confining potential. 



C. The instability of the dressed atom versus the 
transverse field corrections 

By contrast to what happens for Bopt, the imaginary 
part of the scalar potential V(r) does not cancel out. 
This leads to an instability of the "dressed" atom asso- 
ciated with the spontaneous decay of the mixed excited 
state [H, US]. However, the long lifetime Tp^ = 850 ns 
of the ^Pi atomic state contributes to slow down 
the decay of the dressed atom governed by the rate 
F' = {'ih'^n'^/Wlf)Tp for the chosen detuning (see Ap- 
pendix). There is a clear relation between the magni- 
tudes of Bopt and F' since they are both proportional to 
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FIG. 4: The geometric field B3(r) (in light blue) and the scalar confinement potential (in dark red) for the triangular optical 
lattice within the rotated frame. Left : The holes and the bumps of B3(r) form an hexagonal lattice. They can be organized in 
rows alined along the x-axis. Right : This can be used to introduce in H an extra y-dependent W(r) potential leading to the 
total scalar potential Vc(r) acting on the state Iz = h/2 (in jiTlBq" units). For sake of visibility Vc(r) is shifted up by 1 unit. 



fi^, hence to the laser intensity Following the current 
use, the laser intensity is defined in terms of the saturated 
intensity Ig at which il^ = For the Yb transition 

considered, Is amounts to 0.14 mW/cm^ [22| leading to 
a decay rate F' = 1 s~^, obtained for a laser intensity 
/;= 73 mW/cm^ which gives rise to Bopt = 4.824 G. 
Obviously, the laser intensity will have to satisfy a com- 
promise if one wants to keep F' oc J; small as well as the 
parameter r/ cx If'^ which scales the importance of the 
transverse gauge field corrections. Below we show that 
a realistic compromise does exist, allowing one to antici- 
pate a magnitude of corrections at the percent level. 

To get a more precise evaluation of the transverse 
gauge field corrections than the one given in equation 
(P^ . we have calculated numerically the average of the 
two quantities Af — ' ^'^'^ using 

their exact expressions over the relevant plaquette. In the 
rhs of Eq. ([2T|). the second term, once averaged, is found 
to be ten times smaller than the first one, whose average 
leads to Afav = 1.53445 k^, while \B{r)\^^ = T684 Bq. 

From their ratio and the value of t^o = ^ (^^^I^tTAi^) ^ 
12.64 for ^''^Yb and for Bq — 1 gauss. We deduce 
Vav = %A/'a«/l-B(i")lai, = 6.828 B^'^ (for Bq expressed 
in gauss). This gives an improved estimate of r]n,m (Eq. 
((22)) ) in this well specified geometry, assuming that the 
Bloch wave functions moduli are slowly varying over one 
elementary cell. One can reduce r]av to the one percent 
level, by choosing Bq = 24 G. For the total laser intensity 



this leads to // = 365 mW/cm^, shared among the four 
beams, with F' = 0.2 s for the "dressed" ground state 
lifetime. 



D. Extension to a triangular optical lattice 
involving running waves 

To end this section, we would like to present an exten- 
sion of the above analysis to the case a triangular optical 
lattice, generated by a set of three interfering running 
waves, invariant by rotation of 27r/3-multiple angles. The 
triangular geometry has been observed for the first time 
by Grynberg et al. [23] and the properties of cold atoms 
trapped in such a 2D lattice have already been considered 
in [24] but with a spatial configuration different from the 
one adopted here. In contradistinction to what we do, 
the authors have not taken into account the two hyper- 
fine transitions. 

We consider three circularly polarized, running waves, 
which propagate in the x, y plane. It should be noted 
that the previous square lattice built from stationary 
waves cannot be identified with a running wave lattice 
invariant by rotation of 7r/2-multiple angles, since the 
normal refiexion of an optical wave cannot be described 
by a rotation of tt of the wave electric field. Our opti- 
cal triangular lattice is described by the following optical 
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electric field, 



the following y-dependent scalar potential, 



E(r) = £^ e„(C)exp(ik„ - iC(n)), 

k„ — _R(z, 2n7r/3).x, 
e„(0 = ii'(z,2n7r/3).(y + iC5)/V2. 



(31) 



It turns out that the choice C(n) = 2-1111 /i is slightly 
more favorable than C,{n) = 0. To simplify the writing, it 
is convenient to introduce the dimensionless vector field: 
b(r) = — 2-p-E*(r) A E(r). Performing a straightforward 
computation using the above expression of E(r), one ar- 
rives at the rather compact expressions for the b(r) field 
components: 



6, = i ( - cos (VSy - 2C) + cos (|) cos - C) 



6, = |V3sin(|)sin(^-C), 



C)(cos(^) -cos( 



V3y 







(32) 



By contrast to the previous case, bz is independent of the 
beam helicity. Using equation ([55t the optical effective 
magnetic field is given by: 



Btri 
opt 



(33) 



As before, in order to avoid singularities in the geometri- 
cal gauge fields, we have added to b(r) a constant field: 
bo = ^/2x + ^y, leading to a constant magnetic field of 



magnitude B}[ 
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Following the procedure used 



77 hWhf ' 

in Sec. V. B for the square lattice, we have calculated the 
geometric field B3(r) satisfying the formula (fT6|) . On Fig. 
m the Left graphic gives a perspective 3D plot of Bem(r), 
the electromagnetic field equivalent of the geometric field 
B3(r) = ^ • ^ A A3(r), as a function of fcx and ky. The 
geometric phase $geom(l/2), given by the fluxof mB3(r) 
through the plaquette < kx < 47r/3 cancels, but 
through the half-plaquette < kx < 27r/3it amounts 
to $geom (1/2) = 0.384137. Larger values are obtained 
by deforming the plaquette, e.g. for —0.9 < kx < 1.92, 
-0.2 < ky < 2.8, it becomes $geom (1/2) = 0.475858. 
We have calculated the corrections due to the transverse 
gauge fields by following the procedure already described 
for the square lattice. The result rjav — 6.86 {Bq")^'^ 
turns out to be close to the one obtained in this previous 
case (6.83 B^^). 

The geometric magnetic field B^z exhibits regularly 
spaced pairs of positive maxima and negative minima 
forming hexagonal lattices. It appears that the bumps 
and the holes can be organized in rows parallel to the 
X axis. This offers the possibility to bar cold ^^^Yb 
atoms from visiting the connected domains of the 2D- 
space where B^z < 0, by adding to the H Hamiltonian 



W(r) = lyh-fi B*" cos {VSk (y - yo)) . (34) 

The confining potential plotted in red on the right 
graphic of Fig. 2] corresponds to the parameter choice 
11 = 1. and yo — 0. 

Challenging experiments are currently performed by 
overlaying two commensurate triangular optical lattices 
generated by light at the wavelengths of 532nm and 
1064nm 25|. Different 2D lattice geometries have been 
observed by tuning the relative positions of the two lat- 
tices. In this realization the beams are linearly polar- 
ized, so that the periodic optical magnetic field has a 
fixed direction, normal to the lattice, leading to a null 
gauge magnetic field B3(r). Our work paves the way to 
predicting and designing space-dependent optical fields 
which will endow eventually trapped fermions with or- 
bital magnetism, the magnetic spin quantum number 
of the trapped atoms playing the role of their electric 
charge. 



VI. SUMMARY AND PERSPECTIVES 

The purpose of this paper is the study of the geometric 
phases emerging from the non-trivial topology embedded 
in space-dependent vector fields B(r), coupled to the spin 
of a neutral particle. Instead of working directly on the 
geometry of the spin wave functions, we have relied, in 
analogy with our previous work Q , on a "rotated frame" 
approach, where a space-dependent unitary transforma- 
tion W(r) is applied to line up the local B(r) field along 
the spin quantization axis: S ■ ISU = |B| 6*2. Intro- 
ducing the standard Euler angles d{r),(p{r), giving the 
direction of B, the wanted transformation reads: U{r) = 
exp(^ Sz (/)). exp(^ Sy 9). In the "rotated" Hamiltonian 

% = U\t)'HU{t), appears the non-Abelian gauge field 
A{r) = ihW ■ -^U which encapsulates the geometry of 

B(r): n = pV(2M) + V{r) + \B\Sz, the " rotated" 
momentum p being given by the significant expression, 
p = — — A{r). Relying only on the SU{2) Lie alge- 
bra: [5*, S^] = iheijk S'^, A{r) can be written as a lin- 
ear combination of = Sx, — Sy, — Sz, namely 
A{r) = ELi Afc(r). The three vector fields Aj(r) 
involve the gradients of 6{r), (/)(r) (Eq. ([5])), in particular 



the "longitudinal" vector field is A3(r) = cos0(r) 



Or' 



(r). 



Performing the non-Abelian gauge transformation 
^I' — ^ exp(|-Afc(r) S^)'i> leaves the Hamiltonian H invari- 
ant, if A{r) transforms like a standard SU(2) Yang-Mills 
field. However, .4(r) differs from a Yang-Mills field in 
an essential way: it has a vanishing non-Abelian "curva- 
ture" . This means that, unlike SU{2) Yang-Mills fields, 
^(r) is a purely geometric object, not a dynamical self- 
interacting non-Abelian gauge field. 

Pursuing the analogy with Berry's phase, we have con- 
structed a Born-Oppenheimer-like perturbation scheme: 
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% — Ho + Hi, which takes the place of the adia- 
batic approximation with its corrections. The Hamil- 
tonian Hq contains only the "longitudinal" gauge vec- 
tor field A3(r) which is assumed to play the dominant 
role while the "transverse" fields Ai^2(r), incorporated 
in Hi, generate well defined BO corrections, appear- 
ing only to second order. The spin dependence of Ho 
reads Asir)} + js\B{r)\) S,. {At this stage of 

our work, we limit ourselves to a 5 = 1/2 particle in 
a 2D space.) The role of Berry's phase is taken up 
by the Aharonov-Bohm phase <i>/iB(m) given by a loop 
integral of the vector gauge field A3(r) or by the 2D 
flux integral of the associated "geometric magnetic" field: 
B3(r) = sin(0(r)) (^0(r) ^0(r) - ^0(r) £e{r)) z, the 

eigenvalue m of Sz/h standing for the "effective" electric 
charge. Going back to the laboratory frame, the explicit 
evaluation of ^^X^im) is performed in sec. Ill within 

the BO approximation, where H ~ Ho- In a typical ex- 
periment, both interfering paths remain in regions where 
B3(r) = 0. Then, the Feynman phase factor can be fac- 
tored out from the wave function ^'(r) = ^j(r, ® Xi(i") 
and the evolution of the spatial part il>i{Y) ruled out by 
the Hamiltonian ?^o(A3 = 0), is just identical to that 
of "01 (r). Thus, the change of phase is only governed by 
the spin transformation Xi(ri) = lA{vi)xi- At the inter- 
ference point, ri — Y2, unitarity oihl{v) simply implies 
that $51^ (to) — '^ABi'in). In addition, we give an ex- 
plicit method to estimate the transverse field corrections 
to $^b(to). They are given in terms of a set of parame- 
ters rjn,m which can be viewed as the squared amplitudes 
of the first order Hi corrections to the eigenfunctions of 
Ho. If B(r) is a periodic field with a single wave num- 
ber k, we have derived for rin,m the following estimate, 

\rin.,ra\ < 1/8 ((^ik) V (?i7s M| B (r) |„_) ' . 

In order to illustrate the above conceptual develop- 
ments on a concrete physical case, we have applied our 
formalism to cold ^^^Yb atoms dressed by the beams of 
an optical lattice. The dressing optical field involves two 
stationary waves directed along x and y, each beam car- 
rying the angular momentum per photon. The light 
frequency is tuned midway between the two transitions 
connecting the {Qs'^)^So, F = I = \ ground state to 
the hyperfine sublevels F = ^ and | of the {6s6p)^Pi ex- 
cited state. This choice guarantees that the dressed atom 
instability does not affect the effective dc magnetic field 
acting on the atom (see [1^ and Appendix) . A real static 
magnetic field of comparable magnitude is added to avoid 
possible singularities of the optical gauge field. This ge- 
ometry is simple enough to lead to a tractable analyti- 
cal expression of the rcscaled geometric field, fc~^B3(r), 
which depends only on the geometric properties of the 
system Hamiltonian, namely relative directions of propa- 
gation axes and polarizations of the beams but is indepen- 
dent of the light intensity. This latter only determines the 
magnitude of the corrections associated with the trans- 
verse gauge fields via B(r). The explicit expression of 



B3(r) can be casted under the following form: B3(r) = 
— cos(fca:) cos{kx)F(^sm{kx),sm{ky)'jk'^ z, where is a 
positive definite function of x and y. This offers the 
possibility, within the rotated frame, to confine the 
^^^Yb atoms in 2 D domains where i.B3(r) is posi- 
tive by introducing an extra scalar potential yV(r) — 
Tt-Ii Bo fJ, cos(fc x) cos(fc y), as illustrated on Fig. [TJ 

We have added at the end of this section a short pre- 
sentation of the geometric magnetic field B3(r) emerging 
from the geometry of the optical field Bopt generated by a 
triangular optical lattice. We consider a lattice resulting 
from the interference of three circularly polarized running 
waves, invariant by rotation of 27r/3 multiple angles. The 
associated B3(r) field exhibits a 2D hexagonal geometri- 
cal pattern involving bump and hole pairs. Like in the 
case of the square lattice generated by stationary waves, 
it is possible to confine the ^^^Yb atoms outside the B^z 
holes, as shown on Fig. |31 

This work is based on the sole geometric properties 
of the Hamiltonian. Our initial problem, formulated in 
the rotated frame, has been reduced, within a BO-like ap- 
proximation, to that of a spinless particle, having an elec- 
tric charge cx m, interacting with periodic Coulomb and 
magnetic vector potentials, i.e a standard - though non 
trivial - problem of Solid State Theory. What remains to 
be done, with our concrete example, is to fully determine 
the energy bands and the associated Bloch wave func- 
tions and explore directly their geometry. This program 
is beyond the scope of this exploratory paper. 
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Appendix : Optimal detuning of the dressing beam 

When the light detuning is of the order of the hyper- 
fine (hf) splitting, the Hamiltonian which acts on the 
optically dressed ^Sqp^i ^'''-^Yb ground state, receives 
distinct contributions from the two hf lines. We derive 
them by assuming here for simplicity, that the light wave 
is a running wave so that the light angular momentum 
cx = lE* A E/E* • E is spatially uniform. 

There are two remarks which simplify the calculation. 

1) For the 1/2 — )■ 1/2 hf line the electric dipole and the 
spin matrix elements are proportional, (Wigner-Eckart 
theorem). The contribution to the Hamiltonian reads 



2) For the 1/2— !-3/2hf line, one should note that the 
sum over the excited states involving the projection op- 
erator rF=3/2 = J2mp \P = h'^F-){F = h^fl ^dds up 

to 7^F=i/2 to give the unit operator. As a result for very 
large detunings off the two hf lines, the vector coupling 
contributions nearly cancel each other while the scalar 
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ones add up to ^^-j- . Hence the contribution of the sec- 
ond hf hue is ^ £^§^^{1 + ^Q. 

By contrast, for opposite detunings 6i = —S2 — 
Wfif/2h the vector couphng contributions become equal 
and add up together, while the scalar ones subtract. The 
complete Hamiltonian becomes: Hopt — H^^^ + i?'^^ = 

-1^7(1 + 4^^-) - *^r'/2, with B„pt = affe^^' the 

optical field. The imaginary part r'/2 = is 

free of Iz contribution. 

The calculation can be generalized to the case of the 
stationary wave given by Eq. ([^ . When the two detun- 
ings are opposite, the effective Hamiltonian is now given 
by: 



H. 



opt 



8d^ / 



1 



■ E* • E - 



■ I • E* A E 



(35) 



The effective magnetic coupling, associated with the term 
oc I, leads to the following expression for Bopf(r): 



Bopt(i") = - 



8d2 



iE*(r) AE(r). 



(36) 



The calculation of F' is also valid for stationary waves, 
assuming that V(r) c>c E(r)* • E(r) can be replaced by its 
average over one unit cell. 
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